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ABSTRACT 



A new technique for computing the strangeness content of the proton on the 
lattice is described. It is applied to the calculation of the strange quark 
contribution to the proton's spin, specifically to the evaluation of the matrix 
element (ps | sy^y 5 s | ps). Preliminary results are not in disagreement with the 
EMC experiment. 
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The EMC experiment measuring the g l structure function in polarized deep inelastic muon 
proton scattering gave the surprising result that the strange quark contribution to the proton 
spin was substantial 1 . Their measurement, 

g x {x)dx = .126 + .010 + .015 (1) 

^0 

when combined with information from neutron (3 decay and hyperon semileptonic decays 
gives, for the strange quark content of the proton's spin, 

As = -.19 + .07 ( 2 ) 

The error includes both the experimental errors and an estimate of the theoretical uncertainty 
coming from the use of flavor SU(3) symmetry in the analysis. In view of the fact that the 
proton has no strange valence quarks, this is a surprisingly large value when compared to the 
nonstrange quark content inferred from the same measurements 

Au = +.74 + .05 (3) 
Ad = -.51 + .05 



The proton sigma term also has an anomalously large strange quark contribution. An analysis 

,n 2 

(p | ss | p) 



of %N scattering leads to the conclusion 



(p\ uu + dd + ss \ p) 



~ 0.2 (4) 



Both of these results suggest that strange quark operators may have fairly large matrix 
elements between proton states. This violates expectations based on the naive quark model 
and the Zweig rule. This suggests that there are processes for which the valence quark model 
is inadequate, and the effects of the "sea" are substantial. 

It is clearly desirable to carry out a lattice simulation of the relevant matrix elements. 
Lattice calculations are capable of providing first principles theoretical results based on QCD 
alone. If successful, this would provide one of the rare cases in which QCD could explain 
experimental results beyond the ken of the naive quark model. 

The matrix element which expresses the strange quark contribution to the proton's spin, 

Ims^As = (ps\ qsij^s^s I P s ^ ^ 



can be extracted from a lattice evaluation of the three point function 
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(6) 



where ^(x) is a proton field operator constructed from nonstrange quark fields and P + is the 
projection operator on positive energy states. 
The quark line flow of the three point 
function is illustrated in Figure 1. 

The major computational difficulty in 
calculating this three point function lies in 
the fact that the strange quark loop must be 
evaluated at each point on a time slice in 
order that the momentum transfer to the 




^'VMj °P erator ma Y be constrained to Figure 1 Calculation of the strange quark 

be zero. Each such evaluation involves loop correlation with the proton propagator 

solving the Dirac equation on the lattice, 

well known to be a very computationally intensive activity. This is in addition to the 
computer time necessary to create the lattice field configurations and compute the quark 
propagators needed to construct the proton. 

Two observations may make such calculations feasible. One is that it not really necessary to 
compute the closed loop on every site on a time slice. For example, if one were to do such 
computations on alternate sites in each direction (saving a factor of 8 in computation), the 
effect would be to include an admixture of the largest momentum in each direction, n/a, 
along with the zero momentum component. All of such high momentum components are 
suppressed by the exponential decay of the proton propagator. At large Euclidean time, the 
three point function falls off proportional to exp(-E(p)t), where E is the lattice energy of he 
proton state with momentum p, and t is the interval between the proton source point and the 
^"VMj operator. 

The second observation is that the solution to a lattice propagator equation at zero separation 
can be rapidly convergent. Consider the one dimensional Klein-Gordon equation for the 
propagator on a lattice of 2N + 1 points satisfying Neumann boundary conditions: 

<Kn + l) " 2<K«) + <K*-1) " m 2 §(n) = 5(n) 

(7) 

<KA0 = <|>(tf-i) 

(K-AO = H-N + V) 



Its exact solution at the origin is 
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4>(0) 



1 + e 



(1 - 2N)M 



2e M - (2 + m 2 ) + e (X " 2A ^ M (2e M - (2 + m 2 )) 



(8) 



where sinh(M/2) = m/2. This converges exponentially in N, the extent of the lattice, to its 
N — > °° value, 



1 



2e~ M - (2+m 2 ) 



(9) 



This is a general property of the irrelevance of boundary conditions in the large volume limit. 
It can be understood intuitively from the perspective of a hopping parameter expansion, in 
which 2N steps are required to hop from the origin to the boundary and back, giving finite ./V 
corrections proportional to K 2N . 



Each of these observations leads to an approximation which is capable of saving orders of 
magnitude in computation. Furthermore, both of them can be systematically improved. 
Adding additional points to the measurement 
of the operator decreases the contamination 
of the zero momentum component with high 
momentum contributions. Increasing the 
size of the box on which the propagator 
equation is solved improves the accuracy of 
the result exponentially. 



0.17 
0.16 



0.010 1 




0.008 1 




0.006 




0.004 


/ 


0.002 " 




0.000 




-0.002 


f 

/ 


-0.004 


9 




Let us apply these ideas to the evaluation of 
the strange quark fraction of the proton spin. 
To evaluate the Green's functions, we have 
used a set of 16 quenched lattice field 
configurations of size 16 3 x 24 at |3 = 5.7, 
provided at NERSC by the Bernard-Soni 
collaboration. Quark propagators using 
Wilson fermions are available with these 
lattices, so that the principal computational 
task is the calculation of the strange quark 
loop. The strange quark operator is 
evaluated on a time slice 5 units from a 
pointlike proton source, in an attempt to 

balance the need for an isolated proton state against the requirements of low noise. 
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Figure 2 Convergence of the strange quark 
loop at a single site as a function of box size 
/V for several values of the hopping constant K 



An enormous savings in computational time is realized if the strange quark loop is evaluated 
in a box substantially smaller than the entire lattice. However, the value so obtained should 
be an accurate estimate of the value which would be obtained were the entire lattice used. 
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We have explored this issue for six values of K using boxes of size 3,5,7 and 9. The 
conjugate gradient method was used to perform the inversion. The results are shown below 
in Figure 2. For a range of K values from 0.14 to 0.165, the value of the operator q gij^q s 
shows good convergence to a stable value by the time the box size reaches 9. We have 
subsequently used a box with 9 4 sites for each conjugate gradient inversion. The kappa value 
used for the strange quark was K = .162. 



Three values of the hopping constant were 
used for the valence quarks, K = .140, .160, 
and .164. On these lattices, the critical 
value of the hopping constant is K c ~ .170. 
Using a linearized relation between K and 
the bare quark mass 2ma = 1/k - 1/k c , these 
correspond to quark masses of about .1, .2, 
and .6 GeV, respectively. The value of the 
hopping constant used for the strange quark 
in the loop corresponds to a mass of about 
.15 GeV. 

For each of the 16 gauge field 
configurations, the strange quark loop was 
evaluated on a fixed time slice for every 
other site in each direction. This results in 
a mixing of the desired zero momentum 
component of the operator with components 
of momentum zero or nla in each direction. 
The magnitude of this contamination for 
each of the seven high momentum states is 
exp(-(E(p)-E(0)t). A simple free lattice 
propaga tor estimat e of this gives 

exp(-(VM 2 +4d -M) t) where d is 1,2, or 
3, depending on the number of directions 
along which the momentum is nonzero. For 
the present calculation, where M ~ 1.5 and 
t = 5, This gives suppressions of at least 
e 5 - .01. 
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Figure 3 shows the results of the calculation 
for three values of the valence quark 
hopping constant K. The x axis gives the 
location of the time slice on which the 

proton is detected, relative to the time slice on which the strange quark operator is evaluated 
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Figure 3 The strange quark spin fraction of 
the proton As, for three values of the proton 
valence quark hopping constant, 
K= .164, .160, and .140 



- 6 - 

Ideally, the signal would appear when the proton detection slice moves past the strange quark 
operator, and it would rise asymptotically to a stable value. 

From these results it is not possible to infer a quantitative value of the strange quark's 
contribution to the proton spin. However, they do give a strong hint that the strange quark 
spin fraction is negative and substantial. At the lightest value of the valence quark mass, 
corresponding to K = .164, As systematically departs from zero as the time at which the final 
proton is detected grows, becoming 2 
standard deviations below zero at the largest 
values of t. However, no plateau in As has 
been reached. That the departure form zero 
is not purely a statistical fluctuation is 
indicated by the contrast with the negative t 
side of the K= .164 graph in Figure 3, where 
no signal is theoretically expected, and no 
hint of one appears. 

The graph also shows what would be 
required to carry out a quantitative 
evaluation of As. The problem is that below 
t = 5 units from the proton source, which is 
about as far as one can go before the 
statistical errors become prohibitively large, 
the contamination of the proton propagator 
with higher mass states is considerable. 
One can get a quantitative estimate of the 
magnitude of the admixture of high mass 
contributions by directly examining the 
proton propagator. Figure 4 shows the 
proton propagator, together with its single 
proton component estimated by fitting to the 
tail of the curve. At t = 4, about half of the 
propagator's value is due to the 
contributions of high mass contributions. 



In order to extract quantitative information 
from this sort of analysis, one would need 
either proton propagators determined over a 
larger range of time separations, or contribution; the effective mass between each 

propagators constructed from operators pair of time slices 

which project on the single proton state 

much more efficiently. In either case, what would be required would be for the value of the 
correlation determining As to be stable over several time slices. 
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Figure 4 The proton propagator and its 
percentage saturation by the single proton 
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As expected from the naive quark model, there is no indication of any strangeness content in 
a proton whose valence quarks are taken as very heavy. Only for the lightest u and d valence 
quarks in the proton is there some indication of a signal. The appearance and magnitude of 
the signal in the present non-optimized calculation gives reason to believe that the strange 
quark spin content of the proton can be accurately calculated using lattice techniques. 

This calculation could be improved upon in a number of ways. As is usual in lattice 
calculations, the investment of computational resources allows for many improvements, e.g., 
larger volumes and smaller lattice spacing. It is well known that lattices of this size and 
value of (3 are suspect for the study of weak interaction matrix elements. 5 A better 
calculation of As would be performed most sensibly as a part of a larger investigation in 
which other hadronic properties are studied. Obviously, it would be interesting to compare 
results from a simulation of full QCD with results from the quenched approximation. 
There are additional improvements which carry little cost in computation. Improved proton 
sources isolate the proton after fewer time steps than pointlike sources; the use of such 
sources is now standard. As emphasized by El-Khadra et ah, some hadronic quantities, such 
as the •/A|/-T| c splitting, are sensitive to 0(a) terms in the Wilson action. 6 It is very likely 
that the spin content of the proton is also sensitive to these terms. This problem can be 
addressed using the improved cloverleaf action. 7 Finally, we note that a finite 
renormalization of the matrix element is necessary to explicitly compare the lattice result with 
continuum measurements. We have not included this factor here because our results do not 
constitute a quantitative evaluation of As, but it will be necessary in later lattice calculations 
which are to be compared with experiment. 
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